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Abstract 

We include in statistical model calculations the facts that in the nuclear multifragmentation 
process the fragments are produced within a given volume and have a finite size. The corrections 
associated with these constraints affect the partition modes and, as a consequence, other observ- 
ables in the process. In particular, we find that the favored fragmenting modes strongly suppress 
the collective flow energy, leading to much lower values compared to what is obtained from un- 
constrained calculations. This leads, for a given total excitation energy, to a nontrivial correlation 
between the breakup temperature and the collective expansion velocity. In particular we find that, 
under some conditions, the temperature of the fragmenting system may increase as a function of 
this expansion velocity, contrary to what it might be expected. 

PACS numbers: 25.70.Pq, 24.60.-k 
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I. INTRODUCTION 



The determination of the nuclear caloric curve is of great theoretical interest since it may 
help to clarify the physics underlying the breakup of nuclear systems into many fragments, 
i.e., nuclear multifragmentation. For instance, there has been intensive debate on whether 



negative heat capacities should be observed in nuclear systems 



Qflflaaaaa, 



not to 



mention the fundamental question of whether there are any clear signatures of a liquid-gas 
phase transition in nuclear multifragmentation [t], 10, [ll[ Q, Q|- Experimental studies 
have proven to be essential in providing insight into the main properties of the phenomenon 
(see [14j and references therein). The difficulties in extracting the key quantities for the 
problem from experiments have been extensively discussed Ijj]. As a co nsequen ce of these 



difficulties, conflicting experimental observations have been reported 



22 



3, H,y 
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231 ] and, therefore, it has not yet been possible to obtain a clear picture for the nuclear 



multifragmentation process. 

In spite of these uncertainties, many features have been clearly established, such as the 



app earance of an appreciab 
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e co 
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lective radial expansion in central heavy-ion collisions 
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34J . This is intuitively consistent with the results obtained 



by dynamical approaches (see, for instance, 
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371]). in which matter is strongly 



compressed during the first violent stages of the collision and expands afterwards. 

Although statistical models have turned out to be quite successful in explaining many 



properties observed experimentally [5|, [38j, the calculation of this radial flow lies beyond 
the scope of those statistical treatments. Therefore, the radial flow energy is taken as an 
input parameter in these statistical calculations, where it is assumed that its main effect 
is to subtract the energy associated with the radial expansion from the thermal motion 



see, for example, [33|, |39_|). This picture has been criticized by some authors [32J, |40J since 



the non-zero relative velocity between different regions of the system could prevent matter 
within a given region to coalesce at the breakup stage. This effect has been quantitatively 
investigated in ref. [4(3] using the lattice gas model. 
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In this work, we incorporate, in the Statistical Multifragmentation Model (SMM) 
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431 ]. effects associated with the finite size of the fragments in the radial expansion, by 



imposing the constraint that they must lie entirely inside the breakup volume. Although 
the corrections mentioned above and discussed in ref. [40J should also be considered, they 
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will not be addressed here. In sect. HI] we present these modifications to the standard radial 
flow calculations. Their inclusion in the SMM, together with a brief review of this model, is 
performed in sect. IIHI The main results are presented in sect. [TV] and conclusions are drawn 
in sect. [V] 



II. RADIAL COLLECTIVE EXPANSION 



We initially consider fragments as point particles originating from the breakup of a source 
characterized by its mass and atomic numbers A and Z , temperature T, besides its spheri- 
cal breakup volume V = 4ttR 3 /3. If one assumes that matter expands radially with velocity 
u(r), at a distance r from the source's center, the probability that the energy of a fragment 



lies between e and e + de is given by 
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P(e, r)de 



1 



:exp (-[e + S^(r)]/T 



sinh 



T 



ds 



(1) 



where E^ w (r) = ^miu(r) 2 is the radial expansion energy at r, va\ = mAi, where m denotes 
the nucleon mass, and Ai stands for the mass number of the i-th fragment of a partition 
of the system into My pieces. The average kinetic energy of this fragment can be readily 
calculated from the above equation, 



E i (r)= / sP(e,r)de = -T + E j 
Jo * 



W ( r ) 



(2) 



If we now take into consideration that the fragment has a finite size, and that it must lie, 
entirely, inside the breakup volume V, its average kinetic energy may be written as: 



Ei= / E t (r)P c (r)dr , (3) 
Jo 

where Ri stands for the fragment's radius, and P c (r) is the probability that the fragment is 
created at a distance r from the center. 

If we assume that the expansion is irrotational and that the velocity field is given by 



3 



u(r) = 7 



R 



(4) 



where 7 is a constant, and also that the fragments may be formed with equal probability 
at any point inside the sphere of radius R — Ri, P c {r) = 3r 2 /(R — Ri) 3 , the average kinetic 
energy of the fragment is 



2 

flow 



^ Rj 
~ ~R 



(5) 



where f3\ OV! = I7 2 , thus clearly separating the thermal motion and radial expansion contri- 
butions to the kinetic energy of the fragment. Therefore, the total kinetic energy of the Mf 
fragments of the partition is 



^trans = ~(M f - 1)T + £ flow ^ N A>Z A 



A,Z 



Ra 
~R 



(6) 



One should notice that, following ref. 42|, the center of mass motion has been removed 
from the thermal contribution. In the above expression, N^z denotes the multiplicity of 
a fragment with mass and atomic numbers A and Z, and we have defined £fl ow = \ m fi'i ow - 
In the case where the geometric constraints are neglected, so that Ra = in the above 
expression, e flow represents the flow energy per particle, as the sum gives E^ OVf = er ow A . 
One sees that the inclusion of the finite size of the nuclear fragments clearly reduces the 
amount of energy in the radial expansion. In particular, heavy fragments are more affected 
than light ones. Therefore, since it influences the sharing between thermal and collective 
energy in a way that depends on the fragment masses, this correction affects the partition 
modes, and, as a consequence, the values of other physical observables. 

Finally, if we assume that the fragments are formed when the source has expanded to 
(1 + x) of its volume at normal nuclear density, and that the fragments when formed are at 
normal nuclear density, Eq. ([H]) can be rewritten as: 



E, 



where 



E 



trans 



flow 
A,Z 



'-(M f -l)T + J2N A ,zE_ 



flow 
A,Z ) 



£flow^4 



A,Z 



A 



;i + x)a 



1/3 



(7) 



(8) 
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In order to illustrate the magnitude of the corrections, we show, in Fig. [TJ E^/sfi ow 
as a function of the mass number, for Aq = 168, \ = 2, 5, and 9. Comparison with the 
unconstrained results, i.e. E^/ea ow = A, shows that this effect is important, even at very- 
low densities. Therefore, the predictions of the statistical calculations should be modified 
when these constraints are included. 
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Mass number 



FIG. 1: (Color online) E^^/tRow as a function of the mass number, for Aq = 168. 



III. INCLUSION INTO THE STATISTICAL MULTIFRAGMENTATION MODEL 



We briefly recall the main ingredients of the SMM. In it one assumes that the excited 
source under goe s a prompt statistical breakup, subject to strict mass, charge, and energy 



conservation 
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45] 



Aq = Y,N a ,zA, Z = Y,N a ,zZ 

A,Z A,Z 



E%- s - + E* 



AZ 



(9) 



(10) 



Above, Eq S ' represents the ground state energy of the source, E* denotes the total excitation 
energy deposited into the system, and e is the elementary charge. The fragment energies Ea,z 
have contributions from the translational motion, as well as from the nuclear bulk, surface, 
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41 



asymmetry, and Coulomb energies |8j. The latter, is calculated through the Wigner-Seitz 
approximation 



461 ] . More specifically, Ea,z reads: 



E 



A,Z 



3 

-B A<Z + £, AjZ + h AZ +2+ A ' z 



(11) 



We stress that the effects discussed in this work are contained in the changes to the last term 
in the expression above, that were discussed in the previous section. The binding energy of 
the fragments, Ba,z, is calculated using the prescription described in ref. [47]] , whereas the 
remaining terms read: 



E AZ — -C c - 



1 



A 1 / 3 \l + x 



1/3 



E 



A,Z 



e 



A + (V) - - ft) A 2/3 



(12) 



(13) 



and 



(3{T) 



T 2 - T 2 



5/4 



(14) 



T 2 + T 2 _ 

We take for all parameters the same values used in ref. |48j, namely, a Coulomb parameter 
Co = 0.720531 MeV, bulk energy density parameter eo = 16.0 MeV, critical temperature 
T c = 18.0 MeV, and surface energy parameter /3 = 18.0 MeV. One should notice that by 
adding the term associated with the Coulomb energy of the homogeneous sphere in Eq. 
( TTUj) to the Coulomb contributions given by the fragments' binding energies and Eq. ([121) . 
one obtains the Wigner-Seitz expression given in ref. [4lJ. It is also worth mentioning that 
constraints on the center of mass motion are also imposed for each breakup partition, so 
that the total kinetic energy is given by Eq. (J7J). 

The breakup temperature is determined, for each partition, by solving Eq. fflOl) . so that it 
is strongly dependent on the partition mode. As the different terms in the sum appearing in 
Eq. fflUj) are affected in different ways according to the size of the fragments they represent, 
the temperature of the system will change appreciably from the value calculated without 
geometrical constraints. 

The average value of a physical observable O a z is calculated through 



6 



(Oa,z) 



E / °A,z exp 


E{A,Z} / N A,zS A ,Z 


E/ ex P 


Yj{A,Z} f n a,zSa,z 





(15) 



where the sum is performed over all possible partitions {A, Z}f of the nuclear system into 
fragments, and the entropy of fragment (A,Z), Sa,z, is calculated through the standard 
thermodynamic relation 



S = - 



dF 



(16) 



where F is the Helmholtz free energy. Since it depends on the temperature of the frag- 
menting system, the weight of the corresponding mode is also influenced by the constraints 



just described. Ref. 48[ provides a detailed presentation on how empirical information is 
incorporated into F, and we refer the reader to that work for details. Except for the inclu- 
sion of the radial expansion, our SMM calculations follow the description of the Improved 
Statistical Multifragmentation Model (ISMM) presented in that work. 



A. Deexcitation of the primary fragments 

Since most excited fragments are detected after they have undergone secondary decay, 



we have used the Weisskopf treatment described in ref. [49J to estimate these effects on the 
fragment energy spectrum. In this approach, the probability that a compound nucleus, with 
total excitation energy e*, emits a fragment (A,Z), whose mass is Pa,z, is proportional to 

™ re*~!>A,ss-e ( Z]z 

i=o Jo 



where 



s 9a]zVa,zPa,z{ £ ) Pr( £ ~ b A,z - £a]z ~ £ ) „ oX 

In the expression above, bA,z represents the separation energy, g^ z denotes the spin de- 
generacy of the state i, <ja,z is the cross-section of the inverse reaction, e^z stands for the 
excitation energy of the emitted fragment, and p(s*) corresponds to the density of states 
of either the decaying nucleus (CN) or the residual fragment (R). We have used the same 



parameters of ref. 



491 ]. except for the binding energies and the level densities. The former 



are the same used in our SMM calculations, whereas the latter are given by the standard 
Fermi-gas expression (e*)/A = aT 2 , but the excitation energy and the breakup tempera- 
ture are taken as the average values obtained through Eq. (|T5|) for each primordial species. 
Therefore, the density of states 



has a different level density parameter a for distinct primary fragment species. This ensures 
consistency with the population of the excited states in SMM and in the secondary decay 
treatment. 

The final kinetic energy spectrum is generated by a Monte Carlo sample of the possible 
decay channels of the primary distribution. More specifically, the excitation energy of a 
given primordial fragment is selected with probability 



The thermal velocity of the decaying fragment is then assigned according to the Boltzmann 
distribution. 

The radial expansion is incorporated by adding to the velocity a contribution given by Eq. 
(j4]). For consistency, the position of the fragment is uniformly sampled within a spherical 
volume of radius R, which is equal to the breakup volume of the system. We also impose 
the constraint that the fragment must lie entirely inside it. The contribution to the kinetic 
energy due to the Coulomb interaction is estimated by considering the repulsion between the 
fragment and the remaining part of the system. We simply assume that the fragment with 
atomic number Zf is situated inside a sphere of charge (Z — Zf)e, homogeneously distributed 
within its volume. The recoil of this core is taken into account when the corresponding boost 
associated with this binary repulsion is added to the fragment's velocity. 

The selection of a specific channel is made with probability 



where the sum runs over all possible decay channels. We have considered the emission of all 
nuclei from A = 1 to A = 10. 




(19) 




(20) 




(21) 



S 



For the selected deexcitation mode, the relative kinetic energy of the products e < e* — 
&a,z — eJlz * s sampled with weight proportional to Eq. ffT8]) . Their velocities, in the rest 
frame of the decaying fragment, are determined by energy and momentum conservation. 
The excitation energy of the residue is then obtained by energy conservation and it reads 
£ *r = e * ~ ^a,z ~ e az ~ £ - The decay chain is followed until the remnant fragment has a 
negligible amount of excitation energy, i.e. it cannot decay by particle emition. 

This Monte Carlo sample is repeated 10 5 times for each primary species. In the end, the 
multiplicities are weighed proportionally to the multiplicity of the primordial fragments. 




FIG. 2: (Color online) Average flow energy as a function of the average flow velocity. The full 
curve illustrates the results without inclusion of geometrical constraints whereas the dashed one 
corresponds to those with their inclusion. For additional details see the text. 



IV. RESULTS 

To investigate the effects of the collective radial expansion in SMM, we study the A = 168 
Z = 75 system, at a fixed breakup density. In Fig. [21 we show the average flow energy, Eq ow , 
calculated through Eqs. ([S]) and (1TBT) . as a function of the radial velocity /3fl ow (dashed line), 
in a case where the system expanded to three times its volume at normal nuclear density, i. e. 
X = 2. The total available excitation energy of the system was taken to be E* / A = 6 MeV. 
Comparison with the standard unconstrained values, represented in this picture by the full 
line, demonstrates that the inclusion of the geometric constraints dramatically suppresses 
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the amount of energy which may be actually used in the radial collective expansion. One also 
observes that the flow energy reaches a maximum value, of approximately Er ow /A = 1 MeV 
at a value of /3fl ow = /3 max ~ 0.105 close to the maximum possible (when all the energy 
available would appear as radial flow), and then drops to zero as /9fl ow increases further. 
This behaviour may be understood as a consequence of the fact that the total entropy of the 
system diminishes as more and more energy is stored into organized motion, reducing the 
accessible phase space associated with partitions leading to large flow energy values. Since 
the expansion velocity was taken to have a fixed value for all partitions, those which include 
large fragments, and consequently have smaller fragment multiplicities, are clearly favored, 
since they lead to smaller flow energies. 



5.4 5.6 

E / A (MeV) 



5.8 6.0 



FIG. 3: (Color online) Average breakup temperature as a function of the average thermal energy. 
The dash-dotted curve illustrates the results obtained when geometrical constraints are not in- 
cluded. The full and dashed curve are the results when these constraints are included: the full 
(dashed) curve corresponds to velocity flow values below (above) the one leading to a maximum 
flow energy, as depicted in fig. [2j The arrows indicate the direction in which the radial velocity 



increases. 



The changes on the preferred partitions reflect themselves on many observables, such as 
the breakup temperature. Indeed, SMM calculations at fixed breakup volume clearly show 
that the breakup temperature becomes smaller if one simply removes the corresponding 
amount of flow energy from the total excitation energy (see, for instance, 50J and references 



10 




FIG. 4: (Color online) Average primary multiplicity (upper panel) and average mass number (lower 
panel) of fragments as a function of the thermal excitation energy. The curves correspond to the 
same cases illustrated in figj3l 

therein). This is illustrated by the dotted-dashed line in Fig. [3j which displays the breakup 
temperature as a function of the thermal excitation energy, in the case where geometrical 
constraints are disregarded. On the other hand, the constraints associated with the collective 
motion causes the breakup temperature, at a fixed total available excitation energy, to rise 
instead of diminishing, as is also shown in this picture. In this case, the thermal energy is 
defined as the difference between the total available excitation energy and the average flow 
energy, i.e., E t h = E* — E& ovl . This behavior may be explained by the reduction of the 
fragment multiplicity, which leads to larger fragments, moving with less flow energy. The 
requirement of energy conservation, Eq. (fTUl) . then leads to a higher temperature than when 
the constraints are not included. 

The changes on the primary fragment multiplicity and on their average fragment size are 
shown in Fig. H] as a function of the thermal excitation energy. As in the previous plot, 
the dashed-dotted line represents the results obtained without geometrical constraints. The 
inclusion of these constraints causes the fragment multiplicity to drop down as the thermal 
excitation decreases, before the average flow energy reaches its maximum value, i.e., for 
Alow < Anax- Then, for f3^ Qvl > /9 max , it keeps going down while the thermal energy increases 
until it reaches the smallest possible value Mf = 1. As expected, the opposite trend is 
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observed for the average fragment size. 

In spite of the important changes on the observables, the energy spectrum of the particles 
still exhibits a shape which is similar to that expected without constraints associated with 
the finite fragment sizes. Indeed, the circles in Fig. represent the average kinetic energy 
of the primary fragments versus their atomic numbers. The simulation has been carried out 
for E* I 'A = 6.0 MeV, (3f\ ow = ^2e& QW /m, and £fl ow = 2.0 MeV. As shown in that figure, the 
results can be fitted by the linear function Ez = 12.3 + 1AZ MeV. If the energy spectrum 
were interpreted disregarding the geometric constraints, and one wrote E z — §T + £fk,w2Z, 
comparison with the fit would lead to T = 8.2 MeV and £fl ow = 0.7 MeV. However, the 
simulation gives T = 5.9 MeV and, as already mentioned, the expansion velocity corresponds 
to £fl ow = 2.0 MeV. Therefore, the neglect of geometric constraints may lead to important 
uncertainties in the interpretation of the experimental observations. 
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FIG. 5: (Color online) Average frament kinetic energy as a function of the fragments' atomic 
numbers. The lines correspond to a linear fit of the results. 

In order to investigate the influence of the effects associated with the decay of the hot 
primary fragments on the energy spectrum, we have employed the deexcitation treatment 
presented in sect. IIII Al The kinetic energy of the fragments after secondary decay is de- 
picted in Fig. [5] by the triangles. As may be noticed, the slope of the spectrum increases 
appreciably and one may adjust a linear function to reproduce its main trends. Then, one 
finds Ez ~ 16.8 + 2AZ MeV. One sees that the Coulomb repulsion among the fragments 
appreciably affects the slope of the distribution, besides the overall enhancement of the frag- 
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ment's kinetic energy. Nevertheless, the slope is still much smaller than what would be given 
by radial flow alone if the geometric constraints were not taken into account. In spite of the 
great simplifications adopted in our deexcitation treatment, we believe that the main effects 
are included in it, so that more refined decay schemes should not change our conclusions 
significantly. Therefore, our results suggest that the consistent treatment of the geometrical 
constraints are very important in interpreting the experimental observations. 

V. CONCLUDING REMARKS 

We have investigated the effects of the inclusion of geometric constraints due to the finite 
size of fragments in multifragmentation at a fixed breakup volume. Our results show that 
the inclusion of these constraints in SMM lead to qualitative different conclusions on the 
behavior of many physical observables as the system undergoes a radial collective expansion. 
In particular, our results suggest that radial flow alone should not be able to explain a very 
large increase of the fragments kinetic energy as a function of the atomic number. Indeed, 
the simulations presented here show that, for a fixed total excitation energy, the amount 
of energy stored in the radial expansion is strongly suppressed. As a consequence, other 
mechanisms should be considered in order to explain the slopes observed experimentally 
in the energy spectra of the fragments. Thus, we believe that interpretations based on 
statistical calculations in which energy flow is simply removed from the total energy should 
be reviewed. 
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